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INTRODUCTION 


During  the  last  several  years  the  Icing  Test  and  Development  Group  at 
the  Air  Force  Flight  Test  Center  has  been  developing  an  innovative  in-flight 
experimental  technique  for  evaluating  aircraft  icing  problems.  This  technique 
involves  the  production  of  a  cloud  of  supercooled  water  droplets  extending 
several  hundred  feet  behind  a  tanker  aircraft.  The  cloud,  several  feet  in 
diameter,  is  produced  by  injecting  water  through  spray  nozzles  arranged  in 
an  array  suspended  below  and  aft  of  the  tanker  aircraft.  A  test  aircraft 
is  maneuvered  into  an  appropriate  position  behind  the  tanker  so  as  to 
immerse  a  portion  of  the  aircraft  in  the  icing  cloud.  Observations  are  n  ■ 
of  the  buildup  and  shedding  of  ice  on  the  immersed  surfaces. 

Although  icing  clouds  have  been  produced  in  this  fashion  and  builduf 
shedding  of  ice  on  aircraft  surfaces  have  been  observed,  it  is  not  clear  that 
the  cloud  properties  correctly  simulate  the  icing  environment  in  natural  clouds. 

It  is  considered  especially  important  that  the  droplet  size  distribution  in 
the  artificial  cloud  correspond  as  closely  as  possible  to  that  observed  in 
natural  icing  clouds. 

The  droplet  size  distribution  in  the  test  region  of  the  cloud  depends  on  a 
large  number  of  parameters  including  the  details  of  the  spray  nozzle  design,  the 
geometric  arrangement  of  the  array  of  spray  nozzles,  the  altitude  and  speed  of 
the  tanker,  ambient  atmospheric  conditions  and  the  details  of  the  flow  field 
generated  behind  the  spray  nozzle  array. 

Since  a  complete  experimental  simulation  of  all  of  these  factors  in  a 
ground  based  facility  cannot  be  achieved,  a  theoretical  model  which  would  charac¬ 
terize  the  essential  features  of  an  artificial  icing  cloud  would  be  very  useful. 

A  survey  of  the  literature  on  this  subject  has  uncovered  only  a  few  attempts 
to  model  the  development  of  a  droplet  cloud  produced  by  an  array  of  spray  nozzles. 

The  earliest  pertinent  work  appears  to  be  a  model  developed  by  Shapiro,  et  al  (Ref.  1). 


Although  their  analysis  was  strictly  one-dimensional  and  allowed  for  only  unifnrr 
droplet  sizes,  the i r  modeling  of  the  mass,  momentum  and  energy  exchanges  was 
fairly  complete.  Several  year*  later,  Hoffman  developed  a  similar  one-dimensionct 
model  (Ref.  2).  His  model  was  more  general  in  that  it  allowed  consideration 
of  both  liquid  and  solid  particles.  On  the  other  hand,  he  didn't  allow  for 
any  phase  change  and  was  therefore  unable  to  account  for  changes  in  drop  size 
due  to  evaporation.  A  decade  later  Willbank  and  Schulz  (Ref. 3)  published  a 
report  that  included  a  very  detailed  treatment  of  the  mass,  momentum  and  energy 
interchange  and  was  also  applicable  to  nonuniform  droplet  size  distributions. 
However,  their  work  was  also  strictly  one-dimensional .  In  a  more  recent  paper, 
Miller  (Ref. 4)  focused  on  the  question  of  the  freezing  out  of  supercooled  drop¬ 
lets  but  concluded  that  the  required  data  was  unavailable  to  perform  calculations 
Since  no  comparison  between  any  of  these  models  and  experimental  data  has  been 
found,  the  validity  of  these  models  has  yet  to  be  assessed. 

Furthermore,  in  all  of  these  studies  the  flow  fields  were  taken  to  be 
strictly  one-dimensional.  Although  some  important  progress  has  been  made  in 
modeling  the  mass,  momentum  and  energy  interchanges  between  tje  liquid  and  gas 
phases,  no  significant  work  appears  to  have  been  done  on  the  modeling  of  tne 
free  expansion  of  the  cloud  and  its  effect  on  the  droplet  properties.  That 
expansion  has  a  two-fold  effect.  It  importantly  alters  the  mass,  momentum 
and  energy  interchange  betweeen  the  phases  (and  hence  the  droplet  properties  in 
the  test  region).  It  also  provides  a  flow  velocity  normal  to  the  axis  of  the 
cloud  which  (a)  increases  the  cloud  diameter,  and  (b)  provides  a  mechanism  for 
stratifying  the  droplet  size  distribution  in  the  cloud  both  horizontally  and 
radially. 

There  is  clearly  a  lack  of  a  suitable  model  to  describe  the  physical 
characteristics  of  an  icing  cloud  generated  in  flight.  This  report  describes  a 
study  undertaken  to  suggest  an  effective  approach  to  developing  such  a  model. 
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This  study  includes: 

(a)  an  examination  of  the  possibility  of  correlating  existing  experimental 
data  using  an  appropriate  set  of  nondimensional  parameters  (see  Appendix); 

(b)  an  evaluation  of  existing  models  of  mass,  momentum  and  energy  exchange 
between  phases; 

(c)  considerations  regarding  alternative  cloud  modeling  approaches  (including 
the  selection  of  a  preferred  approach); 

(d)  additional  studies  conducted  to  lay  the  groundwork  for  future  work  on 
the  recommended  model . 

Recommendations  for  future  work  are  included  at  the  end  of  the  report. 


ANALYSIS 


A.  General  Considerations 

An  examination  of  photographs,  motion  picture  film  and  available  experimental 
data  indicates  that  the  icing  cloud  consists  of  a  freely  expanding  volume  of 
air,  water  vapor  and  water  droplets  trailing  behind  the  spray  nozzle  array. 

The  primary  features  of  interest  are:  (a)  the  physical  extent  of  the  cloud 
and  the  air  velocity  distribution  within  it,  and  (b)  the  characterization  of 
the  droplet  size,  number,  velocity,  and  temperature  distribution  within  the 
cloud.  In  principle,  these  two  facets  of  the  cloud  are  coupled.  However, 
the  results  of  the  previous  one-dimensional  cloud  model  studies  indicate  that 
the  air  flow  properties  within  the  cloud  are  only  slightly  dependent  on  the 
details  of  the  droplet  histories.  On  the  other  hand,  the  droplet  histories 
are  strongly  dependent  on  the  air  velocity  distribution  within  the  cloud  (see 
especially  Ref.  3).  To  put  this  another  way,  the  droplets  develop  within  the 
aerodynamic  flow  field  generated  by  the  array,  but  the  presence  of  the  droplets 
play  only  a  secondary  role  in  determining  the  overall  character  of  the  aero¬ 
dynamic  flow  field. 

Accordingly,  any  model  should  provide  an  appropriate  description  of  the 
aerodynamic  flow  field  inside  of  which  the  droplets  evolve  in  accordance 
with  prescribed  rules  for  the  exchange  of  mass,  momentum  and  energy  between 
the  liquid  and  gas  phases. 

The  essential  feature  of  the  developing  aerodynamic  flow  field  behind 
a  practical  in-flight  array  is  the  free  turbulent  mixing  between  the  injected 
fluids  and  the  air  flowing  through  and  around  the  array.  In  the  absence  of 
fluid  injection,  the  spray  nozzle  array  produces  a  wake  within  which  the  flow 
velocities  are  smaller  than  the  free  stream  value.  The  conventional  Eulerian 
frame  of  reference  is  adopted  in  this  report  in  which  all  velocities  are 
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measured  with  respect  to  the  array.  In  this  frame  of  reference  the  mean  flow 
field  is  steady.  Since  the  average  velocities  of  liquid  and  air  exiting  the 
spray  nozzles  are  smaller  than  the  flight  velocity,  the  cloud  will  maintain  its 
wake-like  character. 

The  classical  theory  of  wake  flows  deals  primarily  with  the  wakes  behind 
solid  bodies.  Very  simple  and  accurate  equations  are  available  for  the  far 
wake  -  that  is,  the  region  far  behind  the  body  producing  the  wa=ke  (Ref .  5). 

The  theory  for  the  near  wake  (i.e.,  between  the  base  of  the  body  and  the  far 
wake)  is  in  a  much  less  satisfactory  state.  This  is  primarily  due  to  the 
existence  of  large  recirculating  regions  immediately  behind  the  base  of  a 
solid  body.  Such  recirculating  regions  may  be  expected  to  be  quite  small 
behind  the  open  spray  nozzle  arrays  characteristic  of  the  present  problem. 

Although  the  calculation  of  this  type  of  near  wake  should  be  much  easier  than 
for  that  behind  a  solid  body,  no  such  theory  is  available  in  the  literature. 

Any  icing  cloud  model  will  have  to  address  this  problem. 

Before  alternative  models  of  the  icing  cloud  are  described,  a  brief 
digression  regarding  the  far  wake  is  in  order.  The  classical  (i.e.,  no  injection) 
far  wake  analysis  is  based  on  the  observation  that  there  is  a  region  of  decreased 
velocity  in  the  wake  behind  the  body  (Fig.  1).  If  Uq  is  the  axial  speed  of  the 
flow  on  the  axis  and  u  is  the  free  stream  speed,  then  one  can  define  an  axial 
velocity  deficit  by 

Ue  0) 

As  a  result  of  this  velocity  deficit,  there  is  a  corresponding  axial  momen¬ 
tum  deficit  relative  to  the  free  stream  momentum  flux  approaching  the  body. 
According  to  the  principle  of  conservation  of  momentum,  since  there  is  no 
force  on  the  fluid  downstream  of  the  base,  the  total  momentum  defect  in  the 
wake  is  independent  of  the  axial  location,  x,  and  is  equal  to  the  drag  of  the 


body  that  produced  the  wake.  I.e,, 


*_<Jv  =  =  ■jT  ■^6<-p 


(2) 


where  p  is  the  free  stream  density,  R(x)  is  the  local  wake  radius,  r  is  the 
radial  distance  from  the  axis  of  the  wake,  Cg  is  the  drag  coefficient  of  the 
body  and  Ag  is  the  base  area  of  the  body.  For  a  circular  cross  section 

(3) 


where  Rg  is  the  radius  of  the  base.  As  turbulent  mixing  takes  place,  the 
velocity  deficit  on  the  axis  decreases  and  the  wake  radius  increases  such 
that  the  integral  in  eq.  2  remains  constant. 

It  can  be  shown  that  sufficiently  far  downstream  of  the  base  of  the  body 
(the  far  wake)  the  wake  radius  and  velocity  deficit  in  the  wake  are  given  by 


where  p  is  a  turbulent  mixing  coefficient. 

The  same  approach  can  be  used  to  determine  the  asymptotic  (i.e.,  far  wake) 
solution  when  there  is  injeciton.  It  is  instructive  to  do  this  in  two  steps. 

First,  consider  pure  air  injection.  Again,  since  there  is  no  external 
force  acting  on  the  fluid  downstream  of  the  point  of  Injection,  the  total 
momentum  in  the  wake  must  remain  constant.  However,  in  this  case  the  initial 
momentum  is  increased  by  the  momentum  of  the  injected  air.  The  latter  can  be 
shown  to  be  equal  to  the  thrust  produced  by  the  nozzles  injecting  the  air. 
Consequently  eqs.  2,4  and  5  still  apply  provided  the  drag  is  replaced  by  the 
drag  minus  the  thrust  and  Cg  is  replaced  by  an  effective  drag  coefficient 


defined  by 


Dra^. -TTirujt  - -L  A-j  Co  ^ 

A  similar  argument  can  be  made  when  water  injection  Is  added.  However, 

in  that  case  the  water  must  be  treated  as  a  separate  fluid.  The  relative 

velocity  between  air  and  water  droplets  accelerates  the  droplets  to  the  local 

airspeed  some  distance  downstream  of  the  point  of  injection.  According  to 

Newton's  third  law,  the  force  of  the  air  on  the  droplets  is  equal  and  opposite 

to  the  force  of  the  droplets  on  the  air.  According  to  Newton's  second  law, 

this  produces  a  reduction  in  the  momentum  flux  of  the  air.  If  the  rate  of  mass 

« 

addition  due  to  water  injection  is  wi*-  and  the  injection  velocity  is  character¬ 
ized  by  U^- ,  then  by  the  time  the  droplets  reach  the  far  wake,  the  momentum 
transfer  is  *V(  .  Eqs.  4  and  5  still  apply  provided  the  effective  drag 

coefficient  is  further  modified  as  follows 

Doa<^  -  ~ TWust*  +  WW  ^’^,1)  “ 

This  modified  far  wake  analysis  provides  a  very  simple  tool  for  estimating 
the  effects  of  the  array  geometry  and  fluid  injection  parameters  on  the  size 
of  the  cloud  that  is  generated.  For  example,  eq.  4  shows  that  the  cloud  radius 
can  be  increased  by  increasing  the  size  of  the  array  and/or  the  effective  drag 
coefficient.  From  eq.  7  it  can  be  seen  that  the  effective  drag  coefficient  can 
be  increased  by 

a)  using  an  array  with  a  higher  drag  coefficient 

b)  reducing  the  thrust  of  the  injected  fluid.  This  can  be  achieved  by 
reducing  the  air  injection  velocity  but  more  effectively  by  injecting 
in  the  upstream  direction  so  that  the  nozzles  produce  drag  rather 
than  thrust. 

c)  increasing  the  momentum  exchange  between  the  air  and  the  droplets.  This 
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can  be  achieved  by  Injecting  either  laterally  or  upstream  rather  than 
downstream. 

These  results  are  very  useful  for  predicting  the  overall  size  and  shape 
of  the  icing  cloud  but  provide  virtually  no  Information  about  the  droplet 
properties  in  the  cloud.  A  much  more  detailed  model  of  the  cloud  is  required 
for  that  purpose.  Considerations  regarding  the  selection  of  an  appropriate 
modeling  approach  is  provided  in  the  next  section. 

B.  Alternative  Modeling  Approaches 

From  the  previous  discussion  it  is  clear  that  the  modeling  approach  to  be 
chosen  should  be  capable  of  describing  the  development  of  a  wake-like  region 
with  air  and  water  injection  at  the  upstream  end  of  the  wake.  Furthermore, 
it  must  be  capable  of  accounting  for  the  exchange  of  mass,  momentum  and 
energy  between  the  water  droplets,  the  water  vapor  and  the  air  entrained  in 
the  wake.  Fortunately,  the  mass,  momentum  and  energy  exchanges  at  a  point 
depend  only  on  the  local  properties  of  the  flow  field. 

The  modeling  of  those  effects  appear  to  have  been  adequately  treated  in  the 
strictly  one-dimensional  models  referenced  in  an  earlier  section.  A  discussion 
of  that  part  of  the  problem  will  be  presented  in  a  later  section  of  this  report. 

The  primary  new  feature  that  needs  to  be  modeled  is  the  development  of  a 
freely  expanding  cloud  with  internal  phase  change  and  mass,  momentum  and  energy 
exchange.  It  was  pointed  out  in  the  last  section  that  the  freely  expanding 
cloud  is  essentially  a  constant  pressure,  turbulent  wake.  Unfortunately, 
not  enough  is  known  about  turbulent  mixing  to  allow  a  detailed  prediction  of 
the  turbulent  flow  field,  even  without  the  complication  of  liquid  injection. 
However,  some  approximate  techniques  developed  for  other  applications  offer 
the  promise  of  being  able  to  predict  the  mean  flow  property  distributions  in  the 
wake. 
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There  are  three  general  approaches  which  might  be  considered. 

1)  Two-dimensional  finite  difference  calculations.  This  is  the  most 
detailed  and  sophisticated  approach  in  which  the  complete  system  of  partial 
differential  equations  are  integrated  numerically.  These  equations  include  a 
turbulence  transport  equation  which  models  the  local  turbulent  mixing. 

Although  success  has  been  had  with  some  problems  using  this  approach,  its  value 
for  this  application  is  problematical  because: 

a)  it  is  complicated  to  formulate  and  program  and  requires  considerable 
execution  time  on  the  computer. 

b)  the  effects  of  the  water  <.  oplets  on  the  turbulent  transport 
equation  are  unknown.  The  development  of  a  dependable  model 
would  be  a  research  project  in  itself. 

c)  the  uncertainties  in  the  characterization  of  the  properties  of 
the  injected  fluid  entering  the  wake  would  appear  to  make  the 
undertaking  of  such  a  large  computational  task  cost  ineffective. 

2)  Scheduled  mixing  model.  This  is  a  relatively  simple  model  which  has 

been  used  successfully  in  a  number  of  jet  and  wake-like  problems  involving 
chemical  reactions  (Ref.  6).  In  this  approach  the  rate  at  which  the  external 
flow  is  entrained  in  the  mixing  region  is  specified  a  priori.  The  mixing 
region  is  then  replaced  by  a  succession  of  segments  within  which  the 

flow  is  treated  one-dimensional ly.  This  reduces  the  partial  differential 
equations  to  a  set  of  ordinary  differential  equations.  The  initial  conditions 
for  each  segment  are  determined  by  using  the  final  conditions  from  the  previous 
segment  and  assuming  that  the  fluid  at  those  conditions  mixes  instantaneously 
and  uniformly  with  an  appropriate  amount  of  the  surrounding  fluid.  By  con¬ 
straining  this  mixing  process  to  be  in  accordance  with  conservation  of  mass, 
momentum  and  energy,  the  global  properties  of  the  flow  are  correctly  determined. 

This  approach  is  attractive  because  of  its  simplicity  and  because  the  phase 


change  is  not  unlike  a  chemical  reaction.  However,  there  are  two  weaknesses 
to  this  method  which  are  considered  too  important  to  overlook: 

a)  the  shape  of  the  cloud  boundary  is  too  important  to  be  assumed 
a  priori.  It  is  desireable  that  it  be  predicted  as  part  of  the 
calculation. 

b)  the  use  of  the  one-dimensional  approximation  within  each  segment 
may  suffice  for  problems  in  which  only  the  axial  development  of 
the  flow  properties  is  of  primary  importance.  It  does  not  give 
the  radial  variation  of  the  flow  properties.  Although  the 
scheduled  mixing  model  could  be  extended  by  assuming  a  priori 
the  shapes  of  the  radial  distributions  of  the  flow  properties, 
the  approach  doesn't  appear  to  offer  any  advantages  over  the 
next  method. 

3)  Integral  method.  The  integral  method  has  a  long  history  of  successful 
applications  in  viscous  flows.  In  this  method,  the  lateral  variations  in  fluid 
properties  are  approximated  by  simple  analytical  expressions. 

These  analytical  representations  are  constructed  so  as  to  satisfy  the  boundary 
conditions  and  contain  a  number  of  undetermined  parameters.  These  analytical 
formulas  are  then  effectively  substituted  into  the  partial  differential 
equations  governing  the  problem  and  those  equations  are  analytically  integrated 
in  the  lateral  direction.  The  result  is  a  set  of  ordinary  differential  equations 
for  the  undetermined  parameters.  The  solution  of  those  equations  yields  the 
spatial  variation  of  all  flow  properties  as  well  as  the  shape  of  the  boundary 
of  the  mixing  zone.  Since  the  ordinary  differential  equations  are  very  much 
easier  to  solve  numerically  than  the  partial  differential  equations,  the 
computation  is  both  rapid  and  cost  effective.  Furthermore,  the  simplicity  of 
the  equations  allows  the  effects  of  phase  change  and  mass,  momentum  and  energy 
exchange  to  be  included  without  undue  compl ication. 
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The  weakest  link  in  the  application  of  the  integral  method  to  turbulent 
viscous  flows  has  been  in  finding  a  simple,  yet  accurate  model  of  the  turbulent 
mixing  process.  The  traditional  mixing  length  model  has  sometimes  been  found 
to  belackingin  accuracy  when  applied  to  flows  involving  solid  boundaries. 
However,  it  has  fared  much  better  when  applied  to  free  mixing  flows  without 
solid  lateral  boundaries  (Ref.  5).  That  is  precisely  the  situation  that  needs 
to  be  modeled  in  the  current  project. 

As  a  result  of  these  considerations,  the  integral  method  is  considered  to 
be  the  most  appropriate  choice  for  this  problem. 

The  remainder  of  this  report  is  devoted  to  preliminary  studies  supporting 
the  feasibility  of  the  integral  method  approach  and  providing  some  additional 
insights  into  the  physical  features  of  icing  clouds. 

C.  Preliminary  Form  of  the  Flow  Model 

The  most  prominent  qualitative  feature  of  a  wake  flow  is  that  the  wake 
starts  off  with  a  nonuniform  velocity  profile  which  eventually  transforms, 
through  the  action  of  turbulent  entrainment,  to  a  uniform  profile  with  u-iA* 
(Fig.  1).  For  the  purposes  of  illustrating  this  approach  and  generating  some 
qualitative  information  about  the  cloud,  the  following  preliminary  form  of 
the  model  is  presented. 

In  Ref.  3  it  was  found  that  the  presence  of  droplets  has  a  very  minor 
influence  on  the  flow  properties  of  the  air.  For  the  present  purpose,  the 
injection  of  water  droplets  will  be  neglected.  This  gives  no  information 
about  the  time  history  of  the  droplets  (although  more  will  be  said  about  that 
later  in  this  report).  However,  it  does  give  a  reasonable  qualitative  picture 
about  the  size  of  the  cloud  and  the  velocity  distribution  therein.  It  will  be 
seen  that  this  allows  certain  qualitative  features  of  the  droplet  behavior  in 
freely  expanding  clouds  to  be  inferred. 


To  start,  consider  a  plane  normal  to  the  flow  (Fig.  1).  The  rate 


of  mass  flow  across  this  plane  is 
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The  turbulent  mixing  at  the  boundary  of  the  cloud  entrains  air  into  the 
cloud.  According  to  the  classical  turbulent  mixing  length  theory,  the  rate 
at  which  mass  flow  is  entrained  into  the  cloud  is 


<A**1  =  C>4  ^  2.7T  £  d  S 


(9) 


where  a  is  the  entrainment  coefficient,  u^(x)  is  the  velocity  on  the  centerline 
of  the  cloud  and  e,  is  the  arc  length  along  the  surface  of  the  cloud.  Since 
the  cloud  grows  very  slowly,  ds  may  be  approximated  by  dx.  More  precisely, 

As  =  A*  ^ i  Jr  '  -  The  latter  form  will  be  included  in  future  analyses 

but  its  inclusion  is  not  warranted  for  the  present  purpose.  Differentiating 
eq.  8  and  setting  it  equal  to  eq.  9,  conservation  of  mass  requires  that 
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Since  there  are  no  external  forces  acting  on  the  wake,  conservation  of  momentum 
requires  that  the  momentum  defect  in  the  wake  remains  constant.  That  is 
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Eqs.  10  and  11  are  exact.  The  integral  method  is  introduced  by  assuming 
u(x,r)  with  an  appropriate  expression  which  is  explicit  in  r  and  contains  some 
unknown  parameters  which  are  functions  of  x.  Since  the  qualitative  variations 
of  u  with  r  has  a  cosine-like  appearance  (Fig.l),  it  is  assumed  in  this 
preliminary  form  of  the  model  that 
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where  Uq  and  R  are  undetermined  functions  of  x.  This  form  is  seen  to  satisfy 
the  boundary  conditions 
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In  order  to  carry  out  the  rest  of  the  analysis  (which  can  be  done  without 
further  approximations),  it  is  convenient  to  introduce  the  following  non- 
dimensional  variables. 
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Substituting  eqs.  12  and  14  into  eqs.  10  and  11,  performing  the  integration  with 
respect  to  r  and  the  differentiation  with  respect  to  x  yields 
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where 
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Equations  15  and  16  can  be  integrated  exactly  to  yield 
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where  Rg  and  Vg  are  the  initial  values  of  R  and  V. 

It  may  be  noted  that  V  is  the  nondimensional  velocity  defect  on  the  center- 
line.  In  the  far  wake  the  velocity  defect  disappears  (i.e.,  V-+0  as  x-**,). 

-t/3 

It  can  be  seen  from  eq.  19  that  in  the  limit  as  X  ■  Since,  in  that 

-u ?  1/3 

limit,  eq.  18  shows  that  R~ V  ~  x  ,  we  see  that  this  solution  is  consistent 
with  the  classical  far  wake  solution  presented  in  eqs.  4  and  5. 

The  initial  values  of  R  and  V  can  be  related  to  the  specified  mass  flow 
at  the  nozzle  exit  plane,  the  drag  of  the  array  and  the  thrust  of  the  spray 
nozzles.  That  is,  when  the  appropriate  integrals  of  the  assumed  velocity  profile 
(eq.  12)  are  evaluated  at  x=0,  they  must  give  the  specified  values  of  mass 
and  momentum  flux. 

Parametric  studies  of  the  cloud  flow  field  were  made  using  these  results. 

The  nondimensional  form  of  the  velocity  profile  is  shown  at  various  stages  of 
decay  in  Fig.  2.  The  variation  of  centerline  axial  velocity  defect  with  axial 
distance  is  shown  in  Fig.  3  for  various  values  of  Vg.  Note  that  Vg=0  implies 
zero  initial  velocity  defect  and  Vg=l  implies  zero  initial  centerline  velocity. 
The  variation  of  cloud  radius  with  axial  distance  is  shown  in  Figs.  4  and  5. 

Several  features  of  the  cloud  may  be  deduced  from  these  figures.  First, 
we  note  that  the  slopes  of  these  curves  for  large  X  are  in  accordance  with  the 
classical  far  wake  results.  Next,  we  may  note  that  since  the  spray  nozzle 
array  currently  in  use  has  a  radius  of  about  two  feet  and  testing  is  done 
at  a  position  in  the  cloud  about  three  hundred  feet  behind  the  array,  x/Rg-150. 
Figs.  3  and  4  show  that  this  region  is  near  the  beginning  of  the  far  wake 
regime  (constant  slope  on  the  log-log  plots),  regardless  of  the  value  of  the 
initial  momentum  defect.  In  addition.  Fig.  4  shows  that  the  cloud  radius  at 
this  position  is  about  2  to  2H  times  the  original  radius.  This  appears  to  be 
consistent  with  in  flight  photographs  of  the  cloud. 

As  an  aside,  it  should  be  noted  that  the  simple  cosine  variation  of  velocity 
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with  r  is  qualitatively  reasonable  but  quantitatively  only  a  first  approximation. 
It  is  not  surprising  then  that  at  larger  values  of  Vq  (=.7),  some  anomalies 
develop  for  small  values  of  x/R0  (i.e.,  x/Rq<1 ) .  An  improved  form  for  the 
velocity  profile  which  will  better  approximate  the  true  velocity  profile  will 
be  discussed  later  in  this  report. 

As  a  point  of  reference,  calculations  performed  for  a  test  condition  of 
altitude  of  18,000  feet,  u  =300  ft/s  and  typical  water  air  and  injection 
rates  give  a  value  of  Vq».5  . 

Although  water  injection  has  not  been  explicitly  included  so  far  one 
important  feature  regarding  the  behavior  of  water  droplets  in  a  freely  expanding 
cloud  can  be  inferred.  It  had  been  indicated  earlier  in  this  report  that  all 
previous  models  were  strictly  one-dimensional.  In  that  case  all  streamlines 
and  droplet  trajectories  are  straight  horizontal  lines.  In  a  freely  expanding 
cloud  however,  the  streamlines  and  droplet  trajectories  are  curved.  The  radial 
component  of  the  air  velocity, nr,  can  be  determined  from  the  continuity  equation 
in  cylindrical  coordinates: 
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Equation  20  can  be  integrated  formally  with  respect  to  r  to  give 
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Substituting  eq.  12  into  eq.  21  and  carrying  out  the  indicated  operations  provides 
an  analytical  expression  for  .  It  is  convenient  to  introduce  a  nondimensional 
radial  velocity  defined  by 
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The  result  is 
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This  nondimensional  radial  component  of  velocity  is  plotted  versus  the  non- 
dimensional  radial  coordinate,^,  in  Fig.  6.  Two  important  observations  are: 

(a)  the  radial  velocity  is  inward  and  peaks  at  r/R~,5;  and  (b)  radial  veloci¬ 
ties  approach  ten  percent  of  the  free  stream  speed  where  the  axial  velocity  de¬ 
ficit  is  large.  Droplets  injected  at  the  beginning  of  the  cloud  will  begin 
their  trajectories  at  the  injection  angle. 

As  they  move  they  tend  to  curve  into  alignment  with  the  local  air  velocity 
vector.  Heavier  particles  will  change  direction  more  slowly  than  lighter  parti¬ 
cles  because  of  their  inertia.  In  addition,  the  drag  on  a  droplet  depends  on 
its  size.  More  precisely,  the  force  on  a  droplet  is 

F  =  t f  ('/-vjVoj-nV  (24) 

where  V  is  the  (dimensional)  local  air  velocity  vector,  is  the  local  droplet 
velocity  vector,  Cod  is  the  droplet  drag  coefficient  and  rd  is  the  droplet  ra¬ 
dius.  For  the  low  Reynolds  numbers  characteristic  of  the  small  droplets,  Stokes 
law  gives 

^  XM  M 
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where  u  is  the  viscosity  coefficient  for  air.  Denoting  the  mass  of  a  droplet 


by  m<j,  and  substituting  eq.  25  into  eq.  24, Newton's  second  law  requires 

ma  4&  *  in r/ck}  (7-V&I  (26) 

dx 

Since  m,j=  (4/3)7t ^w,eq.  26  can  be  rewritten  as  the  two  scalar  equations 

(27> 
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where  ud  and  vd  are  the  horizontal  and  radial  components  of  droplet  velocity, 
respectively  and  pw  is  the  density  of  water. 

The  main  point  of  this  discussion  can  be  made  from  eq.  28.  It  can  be 


seen  that  the  acceleration  of  the  droplets  in  the  radial  direction  is  propor- 
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tional  to  rd  .  This  means  that  a  20  micron  droplet  experiences  6.25  times  the 


acceleration  of  a  50  micron  droplet.  As  a  result  a  20  micron  droplet  injected 


at  J--.S  at  a  slight  positive  angle  to  the  axis  may  quickly  align  itself  with  flow 
and  end  up  closer  to  the  axis.  On  the  other  hand,  a  50  micron  droplet  may  continue 
along  a  diverging  trajectory  well  beyond  the  test  region.  In  other  words,  a  radial 
stratification  of  the  droplets  will  be  produced  with  the  larger  droplets  near  the 
outer  boundary  of  the  cloud  and  the  lighter  droplets  near  the  axis. 

This  behavior  can  be  corrected  by  designing  the  array  to  inject  a  prepon¬ 
derance  of  smaller  droolets  near  the  edge  of  the  cloud  and  larger  ones  near  the 
axis.  This  or  other  possible  solutions  to  the  stratification  problem  would  of 
course  need  a  cloud  model  such  as  the  one  suggested  here  to  evaluate  the  proposed 
sol ution. 

D.  Additional  Considerations  Regarding  the  Proposed  Model 

It  was  pointed  out  earlier  in  this  report  that  the  cosine  velocity  distri¬ 
bution  velocity  profile,  although  qualitatively  correct,  cannot  be  expected  to 
provide  accurate  quantitative  results  for  all  conditions  of  interest.  An  al¬ 
ternative  form  would  be  an  n'th  degree  polynomial  where  n  is  large  enough  to 
represent  any  real  flow  of  interest. 

Given  the  limitations  of  turbulent  mixing  theory  and  the  fact  that  little 
detailed  information  is  available  regarding  the  initial  velocity  profile  in  the 
cloud,  a  very  large  value  of  n  does  not  appear  to  be  warranted.  In  fact,  since 
no  data  has  been  found  for  the  initial  velocity  profile  in  the  cloud  there  ap¬ 
pears  to  be  little  justification  to  use  a  value  of  n  greater  than  that  needed 
to  assure  that  all  of  the  boundary  conditions  are  satisfied  and  that  the  correct 
mass  and  momentum  flux  into  the  cloud  are  accounted  for. 

The  boundary  conditions  which  must  be  satisfied  are 

ix t  r-c  u-iAo(x)  ^  -  o 
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Before  describing  the  mass  and  momentum  fluxes,  some  reflections  on  the 
cosine  velocity  profile  are  in  order.  First  note  that  in  eq.  12,  the  initial 
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value  of  Ug  must  be  specified.  This  can  be  accomplished  by  specifying  the  enter¬ 
ing  mass  or  momentum  flux.  One  of  the  weaknesses  of  the  cosine  velocity  distri¬ 
bution  is  that  these  fluxes  cannot  be  specified  independently.  In  order  to  ac¬ 
complish  that  goal  another  undetermined  parameter  has  to  be  introduced.  A  suit¬ 
able  polynomial  which  allows  this  to  be  done  and  satisfies  all  of  the  boundary 
and  symmetry  conditions  is 

[0-v)  (3D 

Specification  of  the  entering  mass  and  momentum  fluxes  determines  the  constant 
8  and  the  initial  value  of  V.  This  results  in  a  seventh  degree  polynomial.  It 
is  hard  to  imagine  that  a  higher  degree  polynomial  could  be  justified  given  that 
no  detailed  knowledge  is  available  regarding  the  velocity  profiles  at  the  nozzle 
exits. 

Parametric  calculations  have  been  made  to  evaluate  the  kinds  of  velocity 
profiles  available  from  eq.  31.  Typical  results  are  shown  in  Fig.  7. 

I 

The  analysis  described  in  the  last  section  for  the  cosine  velocity  profile 
has  been  repeated  using  eq.  31.  An  exact  solution  has  been  obtained  for  the 
relation  between  cloud  radius  and  velocity  defect. 

(32) 
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The  differential  equation  for  V(x)  has  been  reduced  to: 
ilf  _ 
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where  R(V)  is  known  from  eq.  32  and  F  is  given  by 
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Eq.  34  can  actually  be  integrated  in  terms  of  elliptic  Integrals  but  the  complexity 
of  the  result  is  not  very  useful.  Future  work  would  Involve  the  numerical  inte¬ 
gration  of  the  equation. 

Most  of  the  discussion  thus  far  has  focused  on  the  fluid  mechanical  aspects 
of  cloud  development.  The  ultimate  goal  of  the  cloud  model  is  to  determine  the 
properties  of  the  water  droplets  as  they  flow  back  through  the  cloud. 

As  the  droplets  progress  through  the  cloud  they  (1)  accelerate  due  to  the 
action  of  the  surrounding  air,  (2)  cool  off  due  to  thermal  conduction  and  evap¬ 
oration,  and  (3)  decrease  in  size  due  to  evaporation.  Thus  the  cloud  is  actually 
a  mixture  of  air,  water  vapor  and  liquid  water  droplets.  The  water  droplets  in¬ 
teract  with  the  air/vapor  mixture  through  the  exchange  of  mass,  momentum  and 
energy.  The  rates  of  exchange  of  these  properties  depend  on  the  local  conditions 
surrounding  each  droplet,  the  size,  temperature  and  speed  of  each  droplet  and  the 
exchange  coefficients  for  mass,  momentum  and  energy  between  phases. 

The  previous  one-dimensional  work  has  indicated  that  the  presence  of  the 
water  droplets  has  a  very  small  effect  on  the  properties  of  the  air  flow.  That 
is,  the  droplets  behave  as  though  they  are  in  a  bath  of  air  and  water  vapor  but 
the  bath  is  not  significantly  affected  by  the  presence  of  the  droplets.  Prelim¬ 
inary  calculations  made  under  this  contract  indicates  that  for  typical  flight 
conditions  this  assumption  may  sometimes  be  marginal.  For  the  conditions  con¬ 
sidered  in  the  previous  section,  the  injected  air  adds  about  ten  percent  addi¬ 
tional  momentum  to  the  cloud.  Likewise  the  momentum  transferred  to  the  water 
droplets  from  the  air  is  about  20%  of  the  same  figure.  The  net  momentum  exchange 
is  considered  to  be  just  about  large  enough  to  warrant  inclusion  in  modeling  the 
aerodynamics  of  the  cloud  but  not  so  large  as  to  justify  more  then  a  first  order 
refinement  of  the  proposed  fluid  mechanical  model.  On  the  other  hand,  this  momen¬ 
tum  exchange  is  the  primary  driving  force  for  the  dynamics  of  the  droplets.  The 
same  qualitative  remarks  apply  with  regard  to  the  transfer  of  mass  and  energy. 
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E.  Mass,  Momentum  and  Energy  Exchange 

While  the  overall  fluid  dynamic  description  of  the  cloud  Is  not  expected 
to  be  sensitive  to  the  details  of  the  mass,  momentum  and  energy  exchange  be¬ 
tween  the  phases,  it  is  expected  that  the  droplet  properties  will  depend  critically 
on  these  mechanisms. 

The  rates  of  change  of  mass,  momentum  and  energy  of  a  spherical  drop  of 
liquid  have  been  extensively  studied.  In  applying  these  results  care  must  be 
taken  that  (1)  the  thermodynamic  properties  of  water  near  freezing  are  appropri¬ 
ately  accounted  for,  and  that  (2)  the  transfer  coefficients  are  appropriately 
represented  over  the  Reynolds  number  range  of  interest. 

A  survey  of  the  available  data  has  been  undertaken.  Very  little  data 
has  been  found  on  the  thermodynamic  properties  of  liquid  water  and  its  vapor 
near  freezing  temperatures.  The  data  of  Ref.  7  are  the  most  complete  found  to 
date  for  the  specific  heats  and  heat  of  vaporization.  Refs.  8  and  9  provide  data 
on  vapor  pressure  of  the  liquid  at  normal  and  supercooled  conditions,  respectively. 

Ref.  3  provides  the  most  complete  description  of  the  mass,  momentum  and 
energy  exchange  rates  which  have  been  found  so  far.  These  include 

mass  transfer  4^  =  (36) 

at  l  - 

where  md  is  the  mass  of  a  droplet,  k*  is  the  mass  exchange  coefficient,  Mv 
is  the  molecular  weight  of  the  vapor,  Xv  is  the  mass  fraction  of  the  vapor  and 
Xy$  is  the  vapor  mass  fraction  at  saturation  conditions. 

momentum  transfer  4^  = -L p  1 V -  I  TT ^  Co  (37) 

where  Vd  and  ?  are  the  droplet  and  local  air  velocity  vectors,  respectively  and 

Cn  is  the  droplet  drag  coefficient, 
d  _ 

energy  transfer  ^ CT'la  )  +  (38) 

where  CVw  is  the  specific  heat  at  constant  volume  of  water,  T  and  Tj  are  the 
local  gas  and  droplet  temperatures,  respectively,  h  is  the  heat  transfer  coef¬ 
ficient  and  hv  is  the  latent  heat  of  vaporization  of  water. 
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Suitable  low  Reynolds  number  correlations  of  the  transfer  coefficients  are 


summarized  in  Ref.  5  . 


momentum 


energy 


kr-  0+ ^V7') 

U*  (U- }«/“£<:  ''’I 

Pa  y 


where  u  is  the  viscosity  coefficient,  &  is  the  Schmidt  number,  Re  is  the  Rey¬ 


nolds  number  of  the  flow  relative  to  the  droplet 

Re  --  (4; 

Cp  is  the  specific  heat  at  constant  pressure  and  Pr  is  the  Prandtl  number. 


In  order  to  keep  computation  time  as  low  as  possible,  a  start  has  been 
made  to  consider  possible  simplifications.  For  example,  the  droplet  trajectories 
are  governed  by  eq.  37.  Consider  the  simpler  one-dimensional  problem. 

,v\A  ^.^-L-pCu-ui)UTr^1'  £oa  (43) 

It  is  of  interest  to  estimate  the  characteristic  relaxation  distance  for  a  drop¬ 
let  to  attain  say  95%  of  the  gas  speed.  The  most  conservative  estimate  is  ob¬ 
tained  using  U=Uoo  and  the  Stokes  relation  for  Cd^  (Cod=24/Re).  Since 
md  =  ( V3)irrc eq.  43  becomes 

~  (.U*-U4)  -  (44) 


Noting  that 


-  U .  iik 

It  A  -u 


and  introducing  the  nondimensional  quantities 

?  =  —  -  — - 

1  V3  U 


eq.  44  becomes 


VA  Ke  '  (47) 

Equation  47  can  be  integrated  easily.  Using  the  most  conservative  initial  con¬ 
dition,  1^(0)  =  0,  one  obtains 


(48) 


For  a  50  micron  droplet  in  a  lOOm/s  airstream  at  0°C,  the  droplet  reaches  95%  of 
the  airstream  speed  after  only  3.3m.  Smaller  droplets  and/or  slower  airstreams 
would  require  smaller  distances.  This  is  only  about  3%  of  the  distance  to  the 
test  region.  This  illustration  suggests  that  this,  and  perhaps  the  other  trans¬ 
port  processes,  take  place  on  a  time  or  distance  scale  short  compared  to  the  length 
of  the  cloud.  The  implication  is  that  some  sort  of  quasi-equilibrium  approxima¬ 
tion  may  suffice  to  predict  droplet  properties  at  the  test  location  in  the  cloud. 
Further  work  is  needed  before  that  hypothesis  can  be  accepted. 

Another  question  which  must  be  addressed  is  the  treatment  of  the  droplet 
size  distribution.  It  is  known  that  a  range  of  droplet  sizes  will  be  injected 
into  the  cloud.  Since  the  evaporation  and  acceleration  rates  are  sensitive  to 
droplet  size  the  model  should  provide  a  method  to  include  a  given  initial  droplet 
size  distribution.  Preliminary  consideration  has  been  given  to  two  possible  ap¬ 
proaches. 

1)  The  one-dimensional  analyses  which  have  treated  this  aspect  of  the  problem 
have  replaced  the  actual  droplet  size  distribution  by  a  set  of  discrete  droplet 
sizes  (Fig.  8).  In  those  analyses  the  flow  velocity  is  independent  of  the  radial 
coordinate  and  the  calculation  is  straightforward.  In  the  case  of  a  freely  ex¬ 
panding  cloud  the  problem  is  complicated  by  the  radial  variation  of  flow  speed 

of  the  air.  This  can  be  accounted  for  by  further  subdividing  each  discrete  size 
of  droplets  into  a  discrete  set  of  Initial  radial  zones  in  which  the  droplets 
enter  the  cloud. 

2)  An  alternative  approach  is  to  recognize  that  at  a  point  in  the  cloud,  the 
change  in  a  droplet's  properties  will  depend  on  its  size,  its  velocity  and  the 
ambient  conditions.  One  can  in  principle  define  a  droplet  distribution  function 
f(>”d»  Vd ;  x,  r)  and  a  droplet  response  function  F(rd,  Vd;  ambient  parameters);  the 
latter  will  give  the  rate  of  change  of  rd  and  Vd  as  a  function  of  the  ambient  para¬ 
meters.  Conceptually,  a  convolution  of  these  functions  will  give  the  spatial  vari¬ 
ations  of  droplet  properties.  In  a  practical  sense,  it  is  expected  that  this 


k 


approach  might  be  useful  only  if  the  quasi -equilibrium  approximations  alluded  to 
above  turn  out  to  be  justifiable. 
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CONCLUSIONS 


An  evaluation  has  been  made  of  the  primary  features  of  a  freely  expanding 
icing  cloud  in  order  to  determine  the  essential  ingredients  needed  to  satisfac¬ 
torily  model  the  cloud.  It  has  been  concluded  that  the  required  model  must  ade¬ 
quately  represent  (1)  a  wake-like  fluid  mechanical  flow-field,  and  (2)  the  exchange 
of  mass,  momentum  and  energy  between  the  droplets  and  the  surrounding  fluid. 

None  of  the  previous  models  for  describing  the  properties  of  droplets  in 
a  flowing  gas  appears  to  be  applicable  to  this  problem.  All  previous  models  are 
strictly  one-dimensional  and  hence  are  incapable  of  describing  some  of  the  essen¬ 
tial  features  of  a  freely  expanding  cloud.  These  features  include,  in  addition 
to  the  radial  variation  of  the  fluid  mechanical  and  droplet  properties,  an  in¬ 
duced  stratification  of  droplet  sizes  (larger  droplets  toward  the  periphery  of 
the  cloud,  smaller  droplets  toward  the  center). 

The  results  obtained  from  the  one-dimensional  models  show  that  although 
the  droplet  properties  depend  strongly  on  their  surrounding  environment,  that 
environment  is  not  strongly  affected  by  the  droplets.  There  appears  to  be  no 
reason  to  reach  a  different  conclusion  regarding  that  one-way  coupling  for  a 
freely  expanding  cloud.  It  would  seem  then  that  the  droplet  mass,  momentum  and 
energy  exchange  can  be  adequately  modelled  using  the  same  equations  used  in  the 
more  complete  one-dimensional  models.  However,  the  properties  of  the  fluid 
surrounding  the  droplets  must  be  modelled  by  a  more  detailed  fluid  mechanical 
model  which  can  appropriately  describe  the  radial  variation  of  the  fluid  pro¬ 
perties. 

Several  possible  approaches  were  examined.  The  most  promising  approach 
appears  to  be  a  turbulent  mixing  length/integral  method  approach.  A  preliminary 
version  of  the  fluid  mechanical  part  of  the  model  has  been  developed.  The  over¬ 
all  size  of  the  cloud  predicted  by  that  model  appears  to  be  in  agreement  with 
in-flight  data.  As  a  result  of  these  encouraging  observations,  development  of 
a  refined  version  of  the  fluid  mechanical  part  of  the  model  has  been  initiated. 
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At  the  same  time  studies  have  been  initiated  to  examine  possible  simplifications 
in  modelling  the  mass,  momentum  and  energy  exchange  between  phases. 

Perhaps  the  weakest  link  in  the  model  Is  the  specification  of  the  proper¬ 
ties  of  the  fluids  entering  the  cloud  from  the  array.  At  the  present  time  it  ap¬ 
pears  that  there  is  not  enough  data  to  precisely  define  the  radial  distribution 
of  flow  properties  entering  the  wake.  The  proposed  model  uses  a  seventh  degree 
velocity  profile  which  matches  the  global  mass  and  momentum  input.  This  is  pro¬ 
bably  adequate  for  the  fluid  dynamic  part  of  the  model.  It  is  expected  that  the 
accuracy  of  the  model  will  be  limited  by  the  information  available  on  the  radial 
distribution  of  droplet  properties  at  the  upstream  end  of  the  cloud. 

In  this  regard,  it  may  be  noted  that  the  proposed  model  is  equally  ap¬ 
plicable  to  modeling  the  entire  cloud  generated  by  the  array,  as  well  as  the 
individual  cloudlets  generated  by  each  spray  nozzle.  If  detailed  information 
were  available  on  flow  and  droplet  properties  at  the  exit  plane  of  a  spray  noz¬ 
zle,  the  proposed  model  could  be  used  to  model  the  development  of  each  cloudlet 
until  the  individual  cloudlets  overlapped.  The  output  of  that  calculation  could 
then  be  used  as  the  input  for  the  cloud  generated  by  the  total  array  of  spray 
nozzles. 

Finally,  recommendations  are  made  in  the  Appendix  regarding  parameters  to 
be  included  in  experimental  measurements  which  would  facilitate  the  correlation 
of  experimental  data. 
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RECOMMENDATIONS  FOR  FUTURE  WORK 


It  is  recommended  that  work  be  continued  on  the  development  of  the 
modeling  approach  described  in  this  report.  The  suggested  order  of  tasks  in 
this  work  is: 

1.  Investigate  the  possibility  of  simplifying  the  mass,  momentum  and 
energy  exchange  equations.  Preliminary  work  suggests  that  the  time  scale 
for  change  of  these  properties  may  be  very  small  compared  to  the  transit 
time  of  a  fluid  particle  from  the  array  to  the  test  region.  In  that  case, 
a  quasi-equilibrium  hypothesis  can  be  invoked  which  might  drastically 
simplify  the  equations.  This  assessment  can  be  made  by  performing  detailed 
parametric  droplet  trajectory  calculations  for  the  flow  field  of  the  pre¬ 
liminary  model  described  in  this  report,  including  the  complete  description 
of  mass,  momentum  and  energy  exchange.  The  results  of  such  calculations 
should  confirm  or  deny  the  validity  of  the  quasi -equi 1 ibri urn  hypothesis. 

2.  Derive  the  full  equations  for  the  proposed  model. 

3.  Program,  debug  and  execute  sample  computer  calculations  for  the 
proposed  model 

4.  Perform  parametric  calculations  for: 

a)  a  variety  of  flight  test  conditions  of  interest 

b)  a  variety  of  radial  droplet  property  distributions  at  the  upstream 
end  of  the  cloud  (to  determine  the  sensitivity  of  the  results  to 
uncertainties  in  the  initial  conditions). 

5.  Employ  the  computer  program  as  a  simulation  tool  to  assess  the  effects 
of  changes  in  array  and  spray  nozzle  design  on  the  effectiveness  of  the 
icing  cloud. 

6.  Compare  predicted  performance  with  experimental  results. 
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In  addition,  it  is  recommended  that  further  work  be  done,  either  in  a 
wind  tunnel  or/and  in  flight  to  more  completely  characterize  the  initial  flow 
and  droplet  property  profiles  at  the  entrance  to  the  cloud.  Similar  work 
should  also  be  done  for  individual  spray  nozzles. 


APPENDIX 


Nondimensional  Parameters 


A  useful  way  of  correlating  data  from  different  experments  is  to  introduce 
nondimensional  parameters.  Two  approaches  were  taken  to  define  the  nondimensional 
parameters  relevant  to  this  problem:  (1)  the  classical  Buckingham  Pi  Theorem  and 
(2)  considerations  regarding  the  essential  physical  features  of  the  problem. 

The  classical  Buckingham  Pi  Theorem  (Ref.  10)  provides  an  algorithm  for 
deducing  a  set  of  nondimensional  parameters  from  a  given  set  of  dimensional  para¬ 
meters.  Examination  of  the  available  reports  on  icing  cloud  testing  yields  the 
set  of  dimensional  parameters  listed  in  Table  I.  According  to  the  Buckingham  Pi 
Theorem,  these  14  dimensional  parameters,  being  comprised  of  the  four  fundamental 
dimensions,  mass,  length,  time  and  temperature  can  be  formed  into  10  nondimensional 
groups.  While  this  provides  a  comprehensive  list  of  nondimensional  parameters,  it 
is  likely  that  only  a  few  of  them  are  of  first  order  importance  for  this  study. 
Consequently,  considerations  regarding  the  essential  physical  features  of  the 
problem  were  employed  to  reduce  the  list  of  nondimensional  parameters  to  those 
which  are  considered  to  be  of  dominating  importance  for  the  range  of  conditions 
of  practical  interest. 

The  overall  shape  of  the  cloud  can  be  expected  to  depend  on  the  retarding 
effect  (drag)  the  array,  the  array  diameter  and  the  mass  and  momentum  flows 
of  the  injected  water  and  air.  The  rate  at  which  the  injected  fluid  mixes  with 
the  surrounding  air  depends  on  the  turbulent  entrainment  factor,  oi.  .  Thus  the 
longitudinal  and  radial  flow  velocities  can  be  expressed  in  nondimensional  form 
as 


’'Wftnr 
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In  addition,  the  properties  of  the  droplets  at  a  point  in  the  cloud  may  be 
characterized  by  the  droplet  size  distribution  at  the  exit  of  the  spray  nozzles, 
the  relative  humidity  of  the  atmosphere,  the  temperatures  of  the  ambient  atmos¬ 
phere,  injected  liquid  and  injected  air  and  the  ambient  pressure.  Thus 


Here  Re  is  the  Reynolds  number  based  on  the  initial  characteristic  droplet 
radius,  *'3,  and  the  initial  relative  speed  between  the  droplets  and  the  air- 
stream.  Note  that  the  problem  of  droplet  breakup  in  and  immediately  downstream 
of  the  nozzle  exit,  while  an  important  factor  in  this  problem,  is  beyond  the 
scope  of  this  project. 

It  can  be  seen  from  these  functional  relationships  that  the  character! zation 
of  the  droplet  properties  depends  on  the  twelve  nondimensional  parameters 
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The  entrainment  coefficient  *1  is  not  likely  to  vary  greatly  over  the  range 

>  I 

of  conditions  of  interest.  The  drag  coefficient  may  vary  significantly  from  one 
configuration  to  another  and  must  be  included  whenever  comparing  results  using 
different  arrays.  It  must  be  expected  that  the  other  parameters  would  play  a 
significant  role  in  any  attempt  to  correlate  the  existing  experimental  data. 

If  data  were  available  from  existing  experiments  to  evaluate  these  non- 
dimensional  parameters,  one  would  expect  that  correlations  of  all  such  data 
could  be  obtained.  Unfortunately  none  of  the  available  reports  were  found  to 
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contain  enough  information  to  permit  evaluation  any  more  than  just  a  few 
of  these  parameters.  Hence,  no  sensible  correlation  could  be  attempted.  However, 
it  seems  likely  that  most  of  the  information  needed  to  calculate  these  parameters 
were  available  to  the  investigators  but  were  not  included  in  the  reports.  It  is 
strongly  recormended  that  reports  on  future  experimental  work  include  suffi¬ 
cient  information  to  calculate  these  parameters. 

There  may  understandably  be  some  concern  regarding  the  possibility  of 
obtaining  a  useful  correlation  of  experimental  da  . a  from  such  a  large  set  of 
parameters.  It  should  be  noted  that  the  work  reported  in  other  sections  of 
this  report  has  established  relations  between  certain  of  these  parameters. 

Further  work  will  certainly  uncover  additional  relationships.  These  nondimensional 
parameters  are  in  fact  the  key  to  combining  future  experimental  and  theoretical 
work  into  a  useful  design  tool  for  icing  studies. 
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TABLE  I 


DIMENSIONAL  PARAMETERS 


Symbol 


Description 


Measured  Calculate 


D 

9 

k 

L 

m 

Q 

P 

T 

u 

0- 

V 

P 

o 


heat  capacity 

X 

drag  force 

X 

body  force  per  unit  mass 

X 

thermal  conductivity 

X 

length 

X 

mass  flowrate 

volumetric  flowrate 

X 

pressure 

X 

temperature 

X 

ve 1 oc i ty 

X 

X 

thermal  expansion 

X 

viscosi ty 

X 

dens i ty 

X 

surface  tension 

X 

Note:  Length  may  be  drop  diameter  (d)  or  initial  wake  radius  (R  ). 
Velocities  may  include  freestream  velocity  (U)  or  local 
velocity  components — i.e.  horizontal  (uh)  and  radial  (ur) . 
Other  subscripts  may  include  (  )  for  injected  air  and  (  )w 
for  injected  water  quantities. 


39 


TABLE  I  (cont'd) 


DIMENSIONLESS  GROUPS  RESULTING  FROM  BUCKI NGHAM-P I  ANALYSIS 


Hi 

B 

D  /  p  u2  it  L2 

ze 

CD 

coefficient  of  drag 

n2 

« 

p  u  L  /  p 

m 

Re 

Reynolds  Number 

n3 

- 

P  /  P  u2  L2 

m 

Eu 

Euler  Number 

ru 

= 

u2  /  g  L 

m 

Fr 

Froude  Number 

n5 

S 

v  cp  /  k 

m 

Pr 

Prandtl  Number 

n6 

- 

g  B  p2  d3  AT  /  u2 

K 

Gr 

Grashof  Number 

n7 

» 

u2  /  Cp  AT 

a 

Ec 

Eckert  Number 

lie 

= 

P  u2  /  p 

dynamic  over  static 

Us 

= 

u  v  /  a 

surface  forces 

n  j  o 

= 

Q  P  /  m 

mass  flow  ratio 

DIMENSIONLESS  GROUPS 


RESULTING  FROM  OBSERVATION 


M  *  (p  u)w  /  p  U 

R  -  R  /  Rq 

V  -  uh  /  U 

W  *  ur  /  U 

X  -  x  /  R0 

6  -  d  /  L 

X  =  LWC  /  p 
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Figure  1  -  Wake/ Cloud  Geometry 


Figure  4  -  Variation  of  Cloud  Radius  With  Axial  Distance 
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Figure  5  -  Variation  Of  Cloud  Radius  With  Axial  Distance 
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